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Abstract 

The unavoidable finite time intervals between the sequential operations needed for performing 
practical quantum computing can degrade the performance of quantum computers. During these 
delays, unwanted relative dynamical phases are produced due to the free evolution of the superpo- 
sition wave-function of the qubits. In general, these coherent "errors" modify the desired quantum 
interferences and thus spoil the correct results, compared to the ideal standard quantum computing 
that does not consider the effects of delays between successive unitary operations. Here, we show 
that, in the framework of the quantum phase estimation algorithm, these coherent phase "errors", 
produced by the time delays between sequential operations, can be avoided by setting up the delay 
times to satisfy certain matching conditions. 

PACS number (s): 03.67.Lx 
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I. INTRODUCTION 



Building a prototype quantum information processor has attracted considerable interest 
during the past decade (see, e.g., jlj]). This desired device should be able to simultaneously 
accept many different possible inputs and subsequently evolve them into a corresponding 
quantum mechanical superposition of outputs. The proposed quantum algorithms are usu- 
ally constructed for ideal quantum computers. In reality, any physical realization of such 
a computing process must treat various errors arising from various noise and imperfections 
(see, e.g., [l], Q, [J). Physically, these errors can be distinguished into two different kinds: 
incoherent and coherent errors. The incoherent perturbations, originating from the coupling 
of the quantum computer to an uncontrollable external environment, result in decoherence 
and stochastic errors. Coherent errors usually arise from non-ideal quantum gates which lead 
to a unitary but non-ideal temporal evolution of the quantum algorithm. So far, almost all 
previous works (see, e.g., |4, have been concerned with quantum errors arising from 

the decoherence due to interactions with the external environment and external operational 
imperfections. Here, we will not be concerned with these two types of externally-induced 
errors, but will focus instead on intrinsic ones. The coherent errors we consider here relate 
to the intrinsic dynamical evolution of the qubits between operations. This has not been 
paid much attention until a recent work in |8|, where a kind of dynamical phase error was 
introduced. It is well-known that a practical quantum computing process usually consists 
of a number of sequential quantum unitary operations. These transformations operate on 
superposition states and evolve the quantum register from the initial states (input) into 
the desired final states (output). According to the Schrodinger equation, the superposition 
wave function oscillates fast during the finite-time delay between two sequential operations. 
In general, these oscillations modify the desired quantum interferences and thus spoil the 
correct computational results, expected by the ideal quantum algorithms without any oper- 
ational delay. 

Two different strategies have been proposed to deal with these coherent errors. One 

Q 

is the so-called "avoiding error" approach proposed by Makhlin et al in Ref. [9(. Its key 
idea is to let the Hamiltonian of the bare two-level physical system be zero by properly 
setting up experimental parameters. Thus the system does not evolve during the delays. 
This requirement is restrictive and cannot be easily implemented for some physical setups 
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of quantum computing e.g., for trapped ions. A modified approach to remove this stringent 
condition was proposed by Feng in Ref. Q, where a pan of degenerate qnantnm states 
of a pair of two-level systems are used to encode two logic states of a single qubit. During 
the delay these logical states acquire a common dynamical phase, which is the global phase 
without any physical meaning. Thus the above dynamical error can be avoided efficiently 
However, this modified scheme complicates the process of encoding information. Another 
strategy to this problem was proposed by Berman et al. ^]. They pointed out that the 
unwanted dynamical oscillations can be routinely eliminated by introducing a "natural" 
phase, which can be induced by using a stable continuous reference oscillation for each 
quantum transition in the computing process. However, this scheme only does well for the 
resonant implementations of quantum computation. The additional reference pulses also 
complicate the quantum computing process and may result in other operational errors. 

We show in this paper that, in the framework of the quantum phase estimation algorithm, 
the coherent phase errors, produced by the free evolutions of the superposition wave functions 
of bare two-level systems, can be avoided simply and effectively by setting up the delay time 
intervals appropriately. The proposed matching condition can be considered a sort of strobed 
operation (with strobe frequencies corresponding to each different transition energy). For 
simplicity, we simplify each quantum algorithm to a three-step functional process, namely: 
preparation, evolution, and measurement. All the functional operations in this three-step 
process are assumed to be carried out in an infinitesimally short time duration, and thus 
only the delays between them, instead of the operations themselves, are considered. The 
effects of the environment decoherence and the operational imperfections are neglected in 
the present treatment. 

The paper is organized as follows. In Section II, we present our general approach with 
the phase estimation algorithm. Sec. Ill gives a few special demonstrations and shows 
how to perform quantum order-finding and quantum counting algorithms in the presence of 
operational delays. Finally, we give a short summary and discussion in Sec. IV. 



II. PHASE ESTIMATION ALGORITHM WITH OPERATIONAL DELAYS 



Our discussion begins with the phase estimation algorithm 
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12 1 and its finite-time 



implementation with some delays. The programs for some of the existing other important 
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quantum algorithms, such as quantum factoring and counting ones, can be reformulated in 
terms of this problem. The goal of the phase estimation algorithm is to obtain an n— bit 
estimation of the eigenvalue exp(i0) of a unitary operation U T , 

Ut |0)t = e* \<f>) T , (1) 

if the corresponding eigenvector \4>)t, and the devices that can perform operations Ut, U t , 
U T , • • • , and U T " , are given initially. Two quantum registers are required to perform this 
algorithm. One is the target register, whose quantum state is kept in the eigenstate \4>)t 
of the unitary operator Ut- Another one, with n physical qubits and called the index reg- 
ister, is used to read the corresponding estimation results. The needed number of qubits 
n in the index register depends on the desired accuracy and on the success probability of 
the algorithm. The most direct application of this algorithm is to find eigenvalues and 
eigenvectors of a local Hamiltonian Ht by determining the time-evolution unitary operator 
Ut = exp(—iHTt/h). The phase estimation algorithm can be viewed as a quantum nonde- 
molition measurement, and can also be used to generate eigenstates of the corresponding 
unitary operator Ut |l4| . 

The ideal quantum algorithm usually assumes that the quantum computing process can 
be continuously performed by using a series of sequential operations without any time-delay 
between them. In reality, a delay between two sequential operations always exists, intro- 
ducing errors that need to be corrected. For simplicity, we reduce the phase estimation 
algorithm to a three-step functional process, namely: initialization, global phase shift, and 
measurement. All the functional operations in this three-step process are assumed to be 
carried out exactly, and thus only the delays between them, instead of the operations them- 
selves, are considered. Such a simplified finite-time implementation of the phase estimation 
algorithm is sketched in Fig. 1. For convenience we distinguish the physical qubit and the 
logic qubit in the index register. The physical qubit is just a two-level physical system and 
the logical qubit is the unit of binary information. Unlike the scheme in ^(J , wherein two 
physical qubits are used to encode one logical qubit, in the present work one physical qubit 
is enough to encode one logical qubit. The symbol \dj)k with a = 0, 1, j, k = 0, 1, ...,n — 1 
means that the kth logical qubit is encoded by the jth physical qubit. \aj) is the eigenstate 
of the bare Hamiltonian of the jth physical qubit corresponding to the eigenvalue E a . 

The quantum phase estimation algorithm with operational delays can be divided into 
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FIG. 1: Quantum phase estimation with operational delays. Note that: 1) there is an operational 
delay r- (m = 1,2) between successive quantum operations on the jth physical qubit, and 2) 
the jth logical qubit is changed to the (n — j — l)th one after the Hadamard gate H and inverse 
QFT F^ 1 . Here, is time delay between the H and L/j," 3 1 operations, while rj 2 -* is the delay 
between {7f,™ 3 1 and F^ 1 . 

three distinct functional steps: 



A. Initialization 



First, we initialize the index register with n physical qubits in an equal-weight superposi- 
tion of all logical states. This can be performed by applying the Hadamard transform to its 
ground state |0)/ = rij=n-i Note that the target register holds an eigenstate \4>)t of Ut 
with eigenvalue exp(i0). Hereafter, the sub-index I will denote the index state, while the 
subindex T refers to the target state. The computational initial state of the whole system is 

f 1 1 2 " _1 1 

l*(0)> = II % \0 s )j ® |0)t = \ k )i ® l^>r, H 3 = — 

[j=n-l ) j ^ Z k=0 ^ Z 

where \k)i = |ao)n-i <8> • • • <8> |on-i)o are ^e number states of the index register, and Hj is 
the Hadamard transform applied to the jth logical qubit. For convenience, in this paper the 
jth logical qubit is changed into the (n — 1 — j)th logical qubit when applying either the 
Hadamard or the (inverse) quantum Fourier transform (QFT). Of course, the order of the 
physical qubits is not changed. 
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After a finite time delay Tj for the jth physical qubit, the initial state |^(0)) of the 
whole system evolves to 

Mr^}) = { fl 4 (f iE ^\Ui)n-i-i + e- iE ^\l 3 ) n ^M ® |0) r , (3) 

Kj=n-1 V ) j 

with Ej and Ej being the eigenvalues of the Hamiltonian for the jth bare physical qubit 
corresponding to the eigenvectors |0j) and respectively. 

B. Global phase shift 

Second, we shift the "global" phase in the eigenvector of the operator U into a measurable 
relative phase. This can be achieved by using the "phase kick-back" technique Q|. Indeed, 
after applying a controlled- U T J operation c — Uj, defined by 

c-U 3 = 11)^(11®^ + 10)^(01 ®/ T , (4) 

to the jth logical qubit, the state |${rj 1 ' ) }) is evolved to 

i^of}) = n (c-uj)^}) 

j=n-l 

= (e-^°^ ) |Oo)n-i + e-^^llo)^) ® 

• • • ® ^ (e- iE ™-^-|0n- 1 )o + e-^-^e^ll^Oo) ® \<f>) T . (5) 

Here and |0)jj(0| are the projectors of the jth logical qubit. It is the identity or 

unity operation. The controlled- U T J operator means that, if the jth logical qubit in the 
index register is in the state \l)j, the 2 J -fold iteration of Ut is applied to the target register. 
The "global" phase in the eigenvector of the operator U T J is changed as the measurable 
relative phases in the states of the index qubits. 

Before the next step in the operation of the algorithm there is another finite-time delay 

(2) 

Tj for the jth physical qubit. During this time interval each physical qubit of the index 
register evolves again freely according to the Schrodinger equation, while the target register 
is assumed to be still in the state \4>)t- As a consequence, the state of the whole system 
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becomes 

|${r,}) = -L(e-^|0 )„_ 1 + e-^'e^ 1 ^|l )„- 1 )® 

• • • <g> (e~ iE °-^ |0„_!>o + e-^— e^^l^^o) ® |0)t, (6) 

with 

r, = rf + rf (7) 

being the total delay before and after the controlled- C/f.™ 3 1 operation. Note that the dy- 
namical phases of the index qubit can be add up before and after this, as controlled- f/f™ 3 
operator is diagonal in the basis of the (n — j — l)th logical qubit of the index register. 



C. Measurement 

Third, we finally apply the inverse quantum Fourier transform (QFT) on the index register 
to measure the phase in the eigenvector of the unitary operator Ut- The inverse QFT, defined 
by the formula 



2 n -l 



QFT- 1 : \k) — F-'lk) = JL £ exp ( -2tt^ ) 



(8) 



can be performed by using the sequential unitary operations F 1 = F^ = 
HqRq^ ■ ■ ■ H n -2 • ■ ■ Ro,n-i ' ' ' Rn-2,n-i-Hn-i, to the corresponding logical qubits. Here, 

/ 1 \ 
10 
1 

\0 e-™' 2k 



I 



j-k,j 



is a two-qubit controlled-phase operation. It implies that the state |l)j of the target jth 
logical qubit will change by a phase exp(— in/2 k ), if the control (j — k)th logical qubit is in 
the state \l)j-k- If the phase <fi can be exactly written as a n-bit binary expansion, i.e., 



= 27T(.0 o ---0 ri _ 1 ) = ^+ 2ri _ 1 



+ ••• + -77"^ 0j = 0, 1, j = 0,1,..., n- 1, 



(9) 



then the expected final output state of the index register, after applying the inverse QFT, 
is the following product state 



= |0n-l)n-l 



10. 



J/J 



|0o)o- 



(10) 



However, the existing dynamical phase error, arising from the free evolution of the physical 
qubits during the delays, may spoil the desired results. For example, measuring the jth 
physical qubit in the computational basis {|0), |1)}, we have 



Ft; 



The expected result is obtained with the following probability 

1 
2 



P 03 = - [1 + cos(A^)] , £^ = E}-Ef}, (12) 



while an error output state |0j © l)j is obtained with the probability P^.@i = 
[1 — cos(AjTj)] /2. Here © refers to addition modulo two. Note that the above proba- 
bility (12) of obtaining the corret result only depends on the total delay time Tj, but not 
directly on the individual time internals t^; m = 1, 2. 
Obviously, if 

r (D = r (2) 

= 0, i.e., for the ideal algorithm realization without any delay, 
one obtains the desired output \4>j)j- While for the realistic case where r^ l \ ^ 0, the 
required quantum inference may be modified, and thus the real output may not be the 
expected one. A worst case scenario is produced if 

AjTj = (2Z + 1)tt, I = 0,1,2,..., (13) 

because the corresponding error-state output is \<pj © which is incorrect. However, if the 
following matching condition 

A, Tj = 2 (I + l)vr, (14) 

is satisfied, one obtains the desired output \4>j)j, and thus the fast oscillation of the superpo- 
sitional wave function is suppressed in the output of the computation. Above, Tj = Tj +Tj 
is the total effective delay time of the jth physical qubit in the algorithm. The condition 
in Eq. (14) is desirable for implementing quantum algorithms with an arbitrary number of 
qubits and includes as a particular case, the less general condition in |8j for the finite-time 
implementation of the 4-qubit Shor's algorithm. 



III. EXAMPLE AND APPLICATIONS 



We now demonstrate the above general approach via a simple example, and show the 
effects of dynamical phases in finite-time implementations of a few quantum algorithms. 
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A. NOT gate eigenvalue 



First, we wish to determine the eigenvalue of the Pauli operator a X} or NOT gate, by 
running the realistic single-qubit phase estimation algorithm discussed above. Assuming 
that the single-qubit target register is prepared into one of the eigenstates 



(15) 



corresponding to the eigenvalues e 1 ^ with <ft = 0, 7T, respectively. According to the above dis- 
cussions, the final state of the index single-qubit register, after the single-qubit measurement 
just performed by Hadamard transform, can be written as 

|*(r)) 7 = i { [1 + e - iAT+l(t '} e- iE ° T |0>/ + [ 1 - e - iAT+i ^} e~ iElr . (16) 




This implies that the probability for the index register to be finally in the state |0)j or |1) 



is 

P (r) = i [1 + cos <fr cos(Ar) + sin <\> sin(Ar)] , (17) 

or 

Pi(r) = - [1 - cos0cos(Ar) + sin0sin(Ar)] . (18) 

If the target register is in the eigenstate of operator a x with eigenvalue +1, i.e., = 0, 
the probability of getting the expected output |0)/ is Pq{t) = 1, if the condition (14) is 
satisfied. However, if the condition (13) is satisfied, the index register will show the error 
output, i.e., |1)/. 



B. Dynamical phase effects in the quantum order-finding algorithm with delays 

Shor's algorithm for factoring a given number A is based on calculating the period of 
the function f(x) = y x mod A for a randomly selected integer y between 1 and A. Once, the 
order r of y mod A is known, factors of A are obtained by calculating the greatest common 
divisor of A and y r ^ 2 ± 1. A finite-time implementation of the order-finding algorithm can 
be translated to the above quantum phase estimation algorithm with delays. Here, the 
unitary operator whose eigenvalue we want to estimate is the unitary transformation U y , 
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with Uy = I, which maps \x) to \yxmod N) and 

U y \uk) = exp y—^-j \ u k), \uk) = ^exp (^y^j \y x modN), k = 0, ...r - 1. (19) 

By the phase estimation algorithm, we can measure the eigenvalue exp(2nik/r) and con- 
sequently get the order r. However, the present target register can not be prepared accu- 
rately in one of the eigenvectors \uk), as the order r is initially unknown. It is noted that 
Y7k=o \ u k)/V^ = |1)> an d |1) is an easy state to prepare. Thus, the algorithm may be run 
by initially generating a superposition of all eigenstates of the operator U y , rather than one 
of them accurately. 

Without loss of generality, we demonstrate our discussion with the simplest meaningful 
instance of Shor's algorithm, i.e., the factorization of N = 15 with y = 7, which had 
been implemented in a recent NMR experiment jig . In this simplest case, the order r is 
the power of two, i.e., r = 2 n , n = 2, and thus the expected phase estimation algorithm 
can measure exactly the n-qubit eigenvalue k/2 n : k = YjjZo kj^i kj = 0, 1. From the 
measurement eigenvalues we can obtain the order r by checking if y r modN = 1. Following 
the corresponding experimental demonstration we need an index register with n = 2 
physical qubits to measure the eigenvalues of the present unitary operator U y , and a target 
register with m = 4 physical qubits to represent the state |1)t = Ylk=o I m A;)t , /2, \uk)r = 
J^. =0 exp (— 2nikx/2 2 ) \7 X mod 15)r/2, which, in fact, is the equal-weight superposition of 
all the eigenvectors of the operator U y : \x)t — > \7 X mod 15)t, x = 0,1, 2, 3, with U y \uk)r = 
exp (2nik/2 2 ) \uk)r- According to the three-step finite-time implementation of the phase 
estimation discussed in the last section, one can easily prove that the whole system is in the 
following entangled state 

l ^ ^ f l / 2jvi2^ ^ k\ \ 1 

l $ {^}) = 2^II{^ [|0;>i-; + exp y—iAjTj H ~~^2 J IWi-ijj^Wkh, (20) 

before the index register is measured by using the inverse QFT. Here, the unimportant 
global dynamical phase factor exp (—2iEjTj) is neglected. 

In the ideal case, i.e., rj = rj 2 ^ = 0, measuring the index register by the inverse QFT 
will, with a probability equal to 1/4, produce the expected output state 

l^out)/ = \ki)i ® \ko) . (21) 
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Simultaneously, the target register will "collapse" into the state of the corresponding ex- 
pected eigenvector \uk). Once a measurement output, i.e., k/2 2 = {2 l k\ + 2°k )/2 2 is known, 
the order is efficiently verified by checking if y l modN = 1 for % = 2 2 /k, 2 ■ 2 2 /k, ...,r. For 
example, if the output is k = 3, i.e., |\l/out)/ = 1 1 1 ) i ® |lo)cb the order can be verified by 
testing y l modN — 1 for i = {2 2 /3, 2 • 2 2 /3, 3 • 2 2 /3 = 4 = r}. Of course, the algorithm fails 
if the output is k = 0, i.e., the target register collapses into the corresponding eigenvector 
\u ). However, these deductions may be modified in a realistic quantum computing process 
where the delays exist, i.e., ^ 0. In fact, one can easily see from Eq. (18) that, after 

applying the inverse QFT, if the target register collapses into the state \uk), the output in 
the index register reads 

o 



i*out)/ = n 



- (1 + e - iA ^+TTik^ + I (1 - e - iA * T i + ™ k i) [l^j 



(22) 



Therefore, the expected state \ki)i Cg) |fco)o is obtained, only if the delays are set up to satisfy 
the matching condition (14). Otherwise, some errors may appear in the index register. 
In particular, an undesirable bit flip error will be produced if Eq. (13) is satisfied. For 
example, if the target register collapses into the state \us)t, the index register generates a 
null 1 0) j- = |0i)i <8> |Oo)o, but not the expected output |3)j = 1 1 1 ) i ® | l )o- 



C. Quantum counting algorithm with operational delays 

Quantum counting is an application of the phase estimation procedure to estimate the 
eigenvalues of the Grover iteration \v\. Il8|. 

G = -Au A^U f . (23) 

Here, A is any operator which maps |0) to J2x=o \ x ) /V^V : , U maps |0) to — 10) and Uf maps 
\x) to (— l)^ x ^\x). This algorithm enables us to estimate the number of solutions to the 
search problem, as the Grover iterate is almost periodic with a period dependent on the 
number of solutions. Indeed, from the following equation 

Gf|*±)=exp(±27rio; l )|*±), I = 0, 1, 2, .., N, (24) 



with |* ± ) = (|Xi) ± i\X )) /y/2, exp(±2mui) = 1 -2l/N±2i^l/N - (l/N) 2 , and {XJ = 
Yjf(x)=i \ x )/^i l-^o) — 2~2f( x )=o \ x ) /^N — I, we see that either u>i or —u>i can be estimated 
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by using the phase estimation algorithm. This gives us an estimation of I, the number of 
solutions. 

In order to explicitly demonstrate how the dynamical phase error reveals in quantum 
counting, we consider the simple case where I = N/4. The expected eigenvalues we want to 
estimate are exp(±7ri/3), corresponding to the target register being kept in the eigenstates 
\^f±). However, in this case the expected out put u\ = 1/6 cannot be expressed exactly 
in a ra-bit expansion. Following Jones et al. Q| and Lee et al. [3], we now adopt the 
ensemble measurement to approximately characterize the final state of the index register. 
The algorithm operates on two registers: a single-qubit index register and the target register 
with m qubits, which are initially prepared in their ground state: \ip{0))i = |0), \iP(0))t = 
Yl°j= m -i |0)r A quantum counting algorithm with delays can also be performed by three 
operational steps: 

1) Applying the Hadamard transform to two registers simultaneously, we have 

= ® \A)t, (25) 

with = #10), = (|0) + |1))/V2, and |^i)t = c + |^ + ) T + c_|^_) r , c± = 

=Fz exp 

2) After the first finite-time delay we apply the controlled operation c— G = /(1|® 
G T + |0)//(0| <g> I T to the state and have 

1^) = Et| [|0)/ + e^- A ^ ) )|l) / ] ® \9j)t, (26) 
j=± V 

After k repetitions of the above operations, the state of the system becomes 

1*3) = J2^= [|0) J + e^ 7r ^- A ^ (1)+r(2)+ -+ r( ^ 1) )|l) J l ® |^) T . (27) 
j= ± V 2 

Above, the controlled operation c — G means that the operation G is applied to the target 
register only when the control qubit is in state |1)/. 

3) After another finite-time delay T^ k \ we apply a second Hadamard transform to the 
control qubit, producing 

k 

|tf 4 > = \ Ys c i [C 1 + e i(2 ^- Ar) ) |0)/ + (1 - e z{2 ^ l - AT) ) ® |^-) r , r = ^ 



j=i m=l 

(28) 
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and then the expectation value of a z is measured to characterize the final state of the index 
register. This corresponds to determining the population difference between |0)//(0| and 
in the state l 1 !^)), and the result can be expressed as 

< frg >i = cos(2nkuJi - At). (29) 

The expected result for the ideal case, i.e., T^ m > = 0, is < a z >j (r) = cos(2irkuJi), and 
the value uji is estimated by varying A; in a manner based on a technique of Kitaev 
For the present problem, if the number of repetitions of the c — G operator is k = 6, 
the measurement result will be expected as < a z >i= 1. This implies that before the 
measurement the control qubit is in state |0) with a high probability. However, in practice, 
operational delays always exist and thus the wave function of the control qubit acquires 
a nontrivial dynamical phase for each delay. As a consequence, the realistic result of the 
measurement is obviously dependent on the total delay time r = Ylm=i We see again 
that the expected result is obtained only if the matching condition (14) is satisfied. 



IV. CONCLUSION AND DISCUSSION 



Ideal quantum algorithms usually assume that quantum computing can be performed 
by continuously applying a sequence of unitary transforms. In reality, when performing 
a practical quantum computations, there are finite time intervals between the sequential 
operations. During these delays, according to the Schrodinger equation, unwanted relative 
dynamical phases are acquired by the superposition wave function of the physical qubit in the 
quantum register. In general, this phase modifies the desired quantum interference required 
for an ideal quantum computer and thus spoils the correct computational results. Note that 
any entanglement between qubits is caused during these delays, and thus resulting coherent 
phase errors can be avoided by simply setting up the total delay times to satisfy certain 
matching conditions. Under these conditions, the relative physical phases in the final state 
of the superposition wave function are deleted. Of course, the dynamical oscillations, due 
to delays, can also be suppressed by trivially setting up individual delays "rj ; m = 1,2..., 
as AjTj = 2mr. The key observation here is that only the total delay time, instead of 
the duration for every delay, needs to be set up accurately to avoid the coherent dynamical 
phase errors. Therefore, only the proper setting up of the total delay is needed for avoiding 
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coherent intrinsic errors. In these implementations, only the free evolution of the physical 
qubits in the index register is considered. 

Co.npa.ed to prev.ous sche.es fl [1 Q for studying siuular problems, our scherue 
presents some advantages. First, it does not require that the Hamiltonian should be equal 
to zero during the quantum register in the idle state (as done in Q] ) . Second, operations 
to force the generation of additional phases to eliminate these phase errors (as done in |sj]) 
are not needed. Finally, our approach does not need to use a pair of degenerate states, 
formed by using two or more physical qubits, to encode a logical qubit (as done in |lO|) for 
transforming the relative phase into a global phase. Therefore, in principle, our proposal 
should allow the implementation of the expected ideal quantum phase estimation algorithm. 

It is worthwhile to emphasize that only the delays between the sequential functional steps 
of quantum computing are considered in the present simplified scheme. The effective dynam- 
ical phases, acquired by superposition wavefunctions of physical qubits during the effective 
delays, may be added up, as the key operation c — Uj in the phase estimation algorithm 
are diagonal in the logical basis of index register. The applied non-diagonal Hadamard gate 
H and inverse QFT operation F' 1 were assumed to be implemented exactly, and thus the 
coherent errors relating to the possible operational delays inside the initialization and mea- 
surements had been neglected. Indeed, the Hadamard gate had been performed exactly by 



using one-step operation [2Cj, and the one-step operational approach had been proposed 21 1 
to exactly implement the QFT. Furthermore, the present scheme for avoiding the coherent 
dynamical phase error is still robust, even if the operational delays inside the initialization 
and measurement are considered. Usually, only a non-diagonal o^-operation is included in 
a three-step process for realizing a Hadamard gate, and in the inverse QFT for measuring 
a physical qubit. Fortunately, it is not required to add up the dynamical phase before and 
after such a non-diagonal <T x -operation in the quantum phase estimation algorithm. In fact, 
the qubit is not in superposition state before (after) the applied cr x -operation in initialization 
(measurement). Therefore, in the framework of the quantum phase estimation algorithm, 
the present strategy for avoiding the coherent phase error is sufficiently robust. This a 
proach can also be used for other quantum algorithms, e.g., Deutsch-Jozsa algorithm 
wherein the key operation is diagonal and the possible non-diagonal operations are also only 
included in the initialization and measurement operations. 



ap- 
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